Composite pulses in NMR as non-adiabatic geometric quantum gates 
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We show that some composite pulses widely employed in nuclear magnetic resonance experiments 
are regarded as non-adiabatic geometric quantum gates with Aharanov-Anandan phases. Thus, we 
reveal the presence of a fundamental issue on quantum mechanics behind a traditional technique. 
To examine the robustness of such composite pulses against fluctuations, we present a simple noise 
model in a two-level system. Then, we find that the composite pulses possesses purely geometrical 
nature even under a certain type of fluctuations. 
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Geometric phases have been attracting a lot of atten- 
tion from the view point of the foundation of quantum 
mechanics and mathematical physics P, H, d, [1]. Re- 
cently, their application to quantum information pro- 
cessing is spotlighted [1, Q , because they are expected 
to be robust against noise. However, the robustness of 
a geometric quantum gate (GQG), which is a quantum 
gate only using geometric phases, is not completely veri- 
fied. Various examinations on this issue have been re- 
ported 0, i, i, 113, ini, d- Blais and Tremblay 
claimed that no advantage of the GQGs exists com- 
pared to the corresponding quantum gates with dynam- 
ical phases, while Zhu and Zanardi Q showed that their 
non-adiabatic GQGs are robust against fluctuations in 
control parameters. 

In this paper, we show that some composite 
pulses widely employed in nuclear magnetic resonance 
(NMR) [l^, |T3| to accomplish reliable operations is re- 
garded as non-adiabatic GQGs based on an Aharonov- 
Anandan (AA) phase [l^, and propose a simple noise 
model in a two-level system. Then, we classify fluctua- 
tions in terms of the robustness of the GQGs. 

An AA phase appears under non-adiabatic cyclic time 
evolution of a quantum system We note that the 

generalization to the non-cyclic case is given in Ref. [3, 
[la |. Let us write the Bloch vector at i (0 < t < 1) 
as n{t){& K^). We denote a state vector given n{t) as 
|n(t))(e C^). Namely, n{t) ^ (n(t)|cr|n(t)), where a = 
*(f7a;, Uy, Uz)- The symbol * means the transposition of a 
vector. Time evolution is described by the Schrodinger 
equation with the Hamiltonian H{t). Note that \n{t)\ = 
1. Hereafter, we denote n(0) as n. We take the natural 
unit system in which h=\. Suppose that \n{l)) = e'^|n) 
(7 e M): n[l) = n. The AA phase 7g is deflned as [1^ 



7g = 7-7d, 



(1) 
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where 



7d 



{n{t)\H{t)\n{t))dt 



(2) 



is a dynamical phase. 

Next, suppose n+ and n are two Bloch vectors sat- 
isfying (a) n+ • n_ — —1 (i.e., (n-|_|n_) = 0) and 
(b) n±{\) — n± (i.e., there exist 7-1- e M such that 
|n±(l)) = e*"*"* |n±). An arbitrary quantum state |n) 
is expressed by \n) ~ a+\n+) -\- a_|n_), where a± — 
{n±\n). We call n± basis Bloch vector corresponding 
to H{t). The initial state \n) is transformed into the fi- 
nal state |n(l)) = a+ e''''+|n+) + a_ e*^- l*^-)- Thus, the 



time evolution operator U aX t 
(t G [0, 1]) is rewritten as 



1 generated by H{t) 



U = e^^+|n+)(n+| + e^^- |n_)(n_|. 



(3) 



Equation ([3]) becomes a quantum gate with a geometric 
phase, when the dynamical component of 7± is vanishing. 
Let us focus on the Hamiltonian for a one-qubit system, 



H{t) = ^uj{t)m{t) ■ a 



(0<t < 1), 



(4) 



which is inspired by a NMR Hamiltonian. In the case 
of NMR, a;(t) and m(t) are the amplitude of and the 
unit vector parallel to a magnetic field, respectively. The 
dynamical phase vanishes when m(t) • n{€) = p7j . 
We note that the integrand in Eq. ^ is rewritten as 
(n(t)|if(t)|n(t)) = (a;(0/4)tr[(m(t) • cT){n{t) ■ ct)\ - 
(cL'(t)/2)m(t) • n(t), where we use tr[iJ(i)] — and 



tr(criCTj) 



2(5. 



This condition has been widely used 



in the experiments on non-adiabatic GQGs [(|. 

A series of pulses, QOajlSO^QO^, has been widely em- 
ployed in the field of NMR for wide band decou- 
pling [3, [13], where /3fe denotes a spin rotation by 
the angle /3 in degree around fc-axis. This is called 
composite pulse and corresponds to the unitary opera- 



tor e 



,/2 



e which is equal to e 



This is generated by the Hamiltonian 

H(t) = 7rm(t) • cr (0 < t < 1), 



»/2 



(5) 
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FIG. 1: Temporal behavior of the basis Bloch vector '(0, 1, 0) 
during the composite pulse 902,180y90i;. (a) without and (b) 
with fluctuations in the control parameters. The fluctuations 
are given by Eq. (fT5)l . where fo = go — 0.1 and ^ = ?7 = 5. 



where 



*(1,0,0) (0<t<l/4) 
m(i)= <; *(0,1,0) (l/4<<<3/4) 
'(1,0,0) (3/4<<<l) 



Hereafter, we will denote to = 0, ti = 1/4, t2 = 3/4, and 
t3 = 1. Various types of composite pulses have been pro- 
posed [H [13, and their usages have been also discussed 
in the context of NMR quantum computing |,18i] . 

Let us examine the time evolution generated by Hamil- 
tonian (O from the view point of non-adiabatic GQGs. 
We choose n± =*(0,±1,0), where ■ n_ = — 1. Then, 
we have the explicit formula 



n±{t)=± 



sm6{t) sm(p{t) 
- sm6{t) cos (t){t) 
cose{t) 



(6) 



where 



0{t) ^ 2TTt - 



7r/2 




{ti<t< ts) 
(otherwise) 



The temporal behavior of n+ on the Bloch sphere is 
shown in Fig.[T]^a). The trajectory n+ is closed. It 



means that |n+(l)) 



o«7+ 



We find that |n±(l)) 



g=F«'r/2|^_j_^ via solving the Schrodinger equation. We 
note that tt is a solid angle surrounded by the trajec- 
tory n^{t). We also find that m{t) ■ n±{t) = at any 
t € [0, 1], and thus the dynamical component is van- 
ishing. Accordingly, we obtain the non-adiabatic GQG, 
U = e-"/2|n+)(n+| + e'''^^\n_){n_\ = e'^'^'y^'^. One 
of the most commonly employed composite pulses turns 
out a non-adiabatic GQG . 

We will classify fluctuations in terms of robustness of 
the composite pulse 90^:180^903;. A noise model will be 
proposed based on a fluctuated closed curve on the Bloch 
sphere. We examine the situation in which the radio- 
frequency (rf) amplitude and phase, and the resonance 
off-set are temporary fluctuated around their aimed val- 



ues. The fluctuated curve is given by 

sin(0(t) -1- /(t)) sin(0(i)+.9(i)) 
n±(i) = ± I - sin(0(i) + /(i)) cos(0(t) + g(t)) 
cos{e{t)+f{t)) 



, (7) 



where we assume that f(t) and g{t) are continuous and 
smooth in [0, 1] [13] and satisfy 



/(to) = 5(io) = 0, f{t3)^g{h)^0. 



(8) 



We will discuss the relevance of f{t) and g{t) to fluctu- 
ations below. The trajectory 'h±{t) is closed under the 
assumption ([5]), as shown in Fig.[T](b). Tnus, we have 



l^^±(l)) 



(9) 



with a phase 7-1- . Generally, 7± includes both the dy- 
namical and the geometric components. We employ this 
noise model in order to ensure the existence of a definite 
AA phase, although we aware of its artificialit;;^ An anal- 
ysis based on a non-cyclic geometric phase [l2| , [l6| may 
be needed for more comprehensive discussions. 

We derive the Hamiltonian generating the time evolu- 
tion corresponding to Eq. ([7]). By differentiating Eq. ([7]) 
with respect to t E (ti-i,ti) (i = 1, 2, 3), we obtain the 
Bloch equation. Then, we find the Hamiltonian in this 
time interval. Hence, the Hamiltonian at < € [0, 1] is 
given by 



H{t) = ^ij{t)m{t) ■ a 



I dgjt) 
2 dt 



(10) 



where 



oj{t) = 2tt + 



We find that 



dm 

dt '' 



cosicj^it) + g{t)) 
rh{t) = I sini^it) + git)) 




m{t) ■ n{t) = 0. 



(11) 



at any t G [0, 1]. The derivative of f{t) is a fluctuation 
of the rf amplitude, while that of g{t) is that of the res- 
onance off-set. A fluctuation of the rf phase is described 
by g{t). From Eq. ([2]), the dynamical component 7d± of 
7± is given by 



1 



7d± 



to 



dgjt) 

dt 



cos[e{t) + f{t)] dt. (12) 



We show that the following two cases exactly lead to 
7d± = 0. Namely, (i) g{t) ^ t) and (ii) f{t) and g{t) 
have a certain symmetric property under time trans- 
lation. The validity of the case (i) is obvious from 
Eq. We focus on the case (ii). We note that 

90xl80j^90:r has several interesting properties under time 
translation: 9{t + 1/2) — 6{t) + tt, for example. We di- 
vide the total time interval /an = {t (z [tQ,t^]} into the 
four intervals, Ii = {t G [^0,^1]}, ^2 = € [ii,l/2]}, 
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FIG. 2: Temporal behavior of the state vector corresponding 
to the basis Bloch vector '(0, 1, 0) during 901,180^90^. The 
initial state vectors are chosen as |n+) = e'''^'^{\0) 
The solid line is the model with the fluctuations. The fluc- 
tuations are described by Eq. (|15p . where /o = (;o = 0.1 and 
^ = 77 = 5. The dashed line is the ideal case, (a) Re(0|n+(t)). 
(b) lm{0\n+{t)). (c) Re(l|n+(t)). (d) Im{l\n+{t)). 



about more realistic control processes [2l|, |22| . Neverthe- 
less, the present discussion is meaningful to understand 
nature of robustness of a geometric phase. 

We directly solve the Schrodinger equation with 
Eq. (fTU|) in order to calculate the geometric component 
of 7±. First, we choose 

f{t) = /o sin[2<Uj(t)], g{t) ^ go sm[2m]u,{t)], (15) 

at t e [U-i, U], where Ui{t) = {t ~ ii_i)/(t, - U-i) and 
^, 77 €E N. The above functions are piecewise smooth in 
[^0,^3] • We show that the temporal evolution of the 
basis Bloch vector *(0,1,0) during the composite pulse 
90xl80y90:r with the fluctuations in Fig.[T];b). This ex- 
ample corresponds to the case (ii) , since Eq. (|14p is sat- 
isfied. We display the temporal behaviors of |n+(i)) and 
\n+{t)) in Fig.m The state vector \n+{t)) is fluctuated 
around \n^{t)), but |n+(i3)) = |n+(t3)). We find that 
7± = T7i'/2. Thus, 7g± 
discuss another example, 



=F7r/2 is confirmed. Let us 



f{t) = fo sin(8<i), g{t) = go sin(87r?7t). 



(16) 



h ^ {t e [1/2, i2]}, and I4 ^ {t e [t2,t3]}. Let us con- 
sider a case when the conditions 

fit +1/2) = fit), |(t + l/2) = g(t), (13) 

are satisfied. The contribution from Ii (12) to 7d± is 
canceled out by that from I3 (14). Thus, this case leads 
to 7d± = 0. Let us consider another case, in which the 
conditions 



are satisfied. We note that /(1/2) = is imposed in 
Eq. (|14p . In this case, the contribution from Ii (12) is 
canceled out by /4 (Is). This cancellation is related to the 
symmetry 6'(1 — t) = —6it) + n. When fit) and git) have 
a certain symmetric property compatible with the pulse 
sequence, the dynamical phase is vanishing. In addition, 
a case (iii) fit) and git) rapidly oscillate with no correla- 
tion, leads to 7d± ~ 0. We can confirm the validity of the 
case (iii) by numerically solving the Schrodinger equation 
with Eq. (jlOp . The case (i) often happens in experiments. 
From Eq. (|10p . one can find fit) is associated only with 
the amplitude of an external controlled field. This quan- 
tity often shows an overshoot or an undershoot before 
settling a desired strength. One can also encounter the 
case (ii) in experiments. A typical example for Eq. (I13p 
may be an oscillating function, as shown in Eq. (jl6|) . A 
linear combination of such oscillating functions leads to 
7(j-|- = 0. Thus, we expect that a lot of rapid oscillating 
fluctuations approximately satisfy Eqs. p^ or p^ . and 
then 7d± ~ 0. The case (iii) is natural when the origins 
of fit) and git) are independent. These three conditions 
lead to 7d± = 0. Thus, the quantum gate under them 
is still regarded as a GQG. It is necessary to examine 



where /o (go) is a positive real number and ^ (77) is an 
integer (to < t < t^). The above functions also satisfy 
Eq. ([5]). Solving the Schrodinger equation numerically 
leads to 7± = 7g± = ^t:/2. The above results mean 
that the solid angle surrounded by ri±(i) is always tt. We 
conjecture that, as long as the fluctuations are introduced 
by Eqs. ([7]) and ([H]), no dynamical phase should exactly 
lead to 7g± = 7g±. 

It is interesting to study the case in which m(t) •n(t) 7^ 
0. Let us consider a simple operation on the Bloch sphere: 
*(0, 0, 1) ^*(1, 0, 0). This process is reahzed by using 



(1^) either e-'"^^ or e^'"^* 



iO<t< 1), where Ha = nay/A 
and Hb = T^io'x + o'z)/2V2. The former satisfies the 
condition m(t) • n(t) = 0, but the latter does not. We 
describe fluctuations in the two models such as Eq. 



HAit) 

HBit) 



TT 

2 ' 

TT 
71 



df_ 
dt 



dt 



rhAit) ■ c 
2 

rriBit) 



dg_a^ 
dt 2 



vr 

71 



dg_ 

dt 



2 



where iriAit) = * (cos(7r/2 -f g(<)), sin(7r/2 + g(f)), 0) 
and rriBit) = * (cosg(t), sing(t), 0). Since /(O) = 
/(I) = g(0) = g(l) = 0, which corresponds to 
Eq. ([S]), the unitary operator generated by HAit) niaps 
*(0,0, 1) *(1,0,0) even in the presence of fit) and 
git). On the other hand, the numerical calculation re- 
veals that the one generated by HBit) maps* (0,0, 1) — > 
*(0.95, -0.26, -0.16) [Fig.[3]. The results mean that 
Eq. ^ does not always ensure robustness in the present 
model. We can flnd an additional term appears in Eq. © 
when mit) ■ nit) ^ 0. Thus, it may cause a strong fluc- 
tuation. We guess that m(t) - nit) = might play an 
important role for stable time evolution in the present 
model. 
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-1 11 1-1 1 

(c) (d) 

FIG. 3: Temporal behavior of the Bloch vector starting from 
'(0,0,1) under the Hamiltonian is shown in (a) and its 
trajectory projected on najTiy-plane is shown in (c). The fi- 
nal point is*(l,0,0). Temporal behavior of the Bloch vector 
starting from '(0,0, 1) under the fluctuating Hamiltonian 
ifo = go = 1-0 and ^ = t; = 10 in Eq. (|16|) ) is show in (b) and 
its trajectory projected on n^^riy-plane is shown in (d). The 
final point is '(0.95, -0.26, -0.16). 



In conclusion, we showed that the composite pulse 
90a;180y902: is regarded as a non-adiabatic GQG. In addi- 
tion, we proposed a simple noise model based on a fluc- 
tuated curve on the Bloch sphere, and then classified 
fluctuations in terms of robustness of 90i:180j,90a;. Al- 
though the present analysis is artificial, it is suitable for 
evaluating errors in non-adiabatic GQGs since a definite 
geometric phase exists even in the presence of fluctua- 
tions. It is important to improve the present method 
in order to examine a more realistic control process or 
a stochastic process. The fluctuations that we discussed 
should be called regular fluctuations, because the fluctua- 
tions are expressed by the two smooth functions f{t) and 
g{t). On the other hand, when fluctuations are given by 
uniform random variables, even a cyclic evolution may 
not be guaranteed [2^ and thus the robustness is not 
expected as discussed in Ref. [7\ . We emphasize that it 
is important to specify fluctuations in order to evaluate 
robustness of a gate. 
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